clc;
syms x y z u(x,y,z);
v=jacobian(u,[x,y,z]);
simplify(curl(v,[x,y,z]))